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Max-min dependence coefficients 
for Multivariate Extreme Value Distributions 

Helena Ferreira Department of Mathematics, University of Beira Interior, Covilha, Portugal 

Abstract: We measure the dependence among sub vectors of a random vector with Multivariate Extreme 
Value distribution by using the expected value of a range and relate this coefficient of dependence with 
the multivariate tail dependence and extremal coefficients. The introduced coefficient extends the concept 
of madogram for several locations and several regions. The results are illustrated with some usual distri- 
ct ' butions and applied to financial data. 
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^ ; 1 Introduction 

in 

CN | The dependence structure of a Multivariate Extreme Value (MEV) distribution is completely char- 

acterised by its dependence function (Resnick (1987), Beirlant et al. (2004)). Since this function 
cannot be easily inferred from data the dependence coefficients are useful, despite the fact that 
one coefficient cannot preserve all the information about this function. 

The most popular of the dependence coefficients are those based on the tail dependence (Sybuya 
(1960), Li (2009)). They resume the probability of occurrence of extreme values for one or more 
random variable given that another(s) assumes extreme values too. For the MEV distributions the 
extremal coefficient (Tiago de Oliveira (1962-63), Smith (1990)) is certainly a crucial and perhaps 
insurmountable tool when we have to summarize the dependence. For a d— dimensional random 
vector we have 2 d — d extremal coefficients which consistency properties are discussed in Schlather 
and Tawn (2002). For an overview of other dependence measures for dependence see, for instance, 
t^J- ! Joe (1997). 

' To the best of our knowledge, there is no extremal dependence coefficients for p > 3 subvectors 

ON . Xi, ...,X p of a random vector X = (Xi, ...,Xd) with MEV distribution. 

The need of to evaluate the strength of dependence among sub vectors arises for instance in the 

setting of max-stable random fields. Let {^i}igR2 be a max-stable random field and 7i, ...I p sets 
■ p 

of locations in R 2 . The joint distribution of Xi, i € Ij, is a MEV distribution and we want to 

3=1 

resume the dependence among the grouped values {X\, i £ Ij} at different regions Ij, j = 1, ...,p. 
This problem is treated by several authors for two variables X\ and X- } corresponding to two 
locations i and j (see Naveau et al. (2009) and references therein) and the obtained results are 
d . extended for two regions I\ and I2 in Fonseca et al. (2012). 

In finance, we are frequently interested in assessing the dependence among several big world mar- 
kets, considering each one as a random subvector. For an application with grouped financial stock 
markets see for instance Ferreira and Ferreira (2012a). 

We propose to evaluate the degree of dependence among sub- vectors Xi, ...,X p of X with MEV 
distribution by using an expected range, which will be referred as a "max-min coefficient". It is 
a summary measure that takes into account the whole group of the extremal coefficients ex of 
Xj, j = 1, ...,p. Our approach is an extension of the modeling for pairwise dependence throughout 
the madogram (Poncet et al. (2006), Naveau et al. (2009)), an extreme- value analogue of the 
variogram (Cressie (1993)), since it enables to resume the spatial dependence structure for several 
locations or regions of locations. 

The proposed moment-based dependence tool takes into account the spread and dependence among 
the subvectors and can be easily estimated. 
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The paper is organized as follows. We introduce in section 2 the dependence coefficient whose is well 
defined for any random vector with MEV distribution, is a function of its copula and is invariant 
with respect permutations of the variables. Its relations with the multivariate tail dependence 
and the extremal coefficients are presented. Based on the expected range coefficient considered we 
compare a MEV distribution with others more concordant distributions and state some bounds. 
In section 3, we compute the max-min coefficients for the marginal distribution of the Multivariate 
Maxima of Moving Maxima process and the Symmetric Logistic distribution. We refer briefly an 
estimator for the max-min coefficients and apply it to grouped financial stock markets. 

2 Max-min dependence coefficients 

Let X = (Xi, Xd) be a vector of unit Frechet random variables, that is, with marginal distribu- 
tion function F(x) = cxp(— x^ 1 ), x > 0, and G denote the Multivariate Extreme Value distribution 
of X. 

The tail dependence function (Huang (1992), Schmidt and Stadtmuller (2006)) of G is defined by 
l(xi, Xd) = — log G(a;^ 1 , x^ 1 ), (x\, Xd) & [0, oo) d . It is a convex function, homogeneous of 
order one and satisfies 

d d 

Y xj < i(xi, ...,x d ) < y^zj, 

with the lower bound corresponding to X with totally dependent margins and the upper bound to 
X with independent margins. 

The tail dependence function l(x\, Xd) gives us information about the probability of occurrence 

d 

of extreme events for the maximum \l F(Xj) given that one fixed margin assumes extreme values 

3=1 

too. In fact, we have (Schmidt and Stadmuller (2006), Ferreira and Ferreira(2011)) 
l(x u ...,x d )=lim(-tlogP(F(X 1 ) <l-^-,...,F(X d ) < 1-y)) = 
lim tP (f(X 1 ) > 1 - ^ V ... V F(X d ) >l-^f) 

and, in particular, 

Id 

l(x, x) = xl(l, 1) = x lim P V F(Xj) > 1 

t— too \ " 

V=i 

For Si(S) = liiieS and 5 t {S) = if i <£ S, it holds that 

l(5 1 (S),...,S d (S)) = e Xs , (1) 

where ex s is the extremal coefficient of the subvector X5 of X with indices in S (Tiago de Oliveira 
(1962-63), Smith (1990)). It takes values in [1, with ex s = 1 when X5 has the minimum copula 
C*x s ••■) u d)s = A u j an d £x s = l^l when X5 has the product copula Cx s (ui, Ud)s = 

jes 

Let X = {Ji, I p } be a partition of D = {1, d}, M(Ij) = \J Xi and X/j the subvector of X 

with indices in Ij. We resume the extremal dependence among X/., j = l,...,p, by the coefficient 
i?(X, A, I) defined as follows. 



i|*W>i-i 
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Definition 2.1. Let X be a vector of unit Frechet random variables and Multivariate Extrem Value 
distribution. For each A = (Ai, A p ) S (0, oo) p and partition I of D, we define 

R(X, A, I) = E I \/ F X > (M^)) - A pX ' ( M &)) 



We remark the relations 



i?(X,A,X) = J B( \/ (^(M^-^Wi))!] (2) 

,{i,j}C{l,...,p} 



and 



i?(X, A,J) =E V \/ F^(^) - /\V J^j(x 4 ) . (3) 

By taking p = 2 in ©, we find in ^i?(X, A, I) the generalized madogram introduced in Fonseca 
et al. (2012), which in turns is the the A-madogram (Naveu et al. (2009)) when d = 2 = p and 
1 - A 2 = Ai e (0,1). 

The max-min coefficient and the generalized madograms for pairs of sets Ii and Ij can be related 
throughout 

i?(X,A,Z)> \/ ^((X^Xi,), (A*, fy). 

{i,j}C{l,-,p} 

d 

We first present a key result that relates the expectation of \J F Xj (Xj ) with the tail dependence 

j=i 

function of G, which enables the derivation of the main properties of i?(X, A, I). The result also 
points out that in this work we can assume that the MEV distributions have unit Frechet margins 
without loss of generality. 

Proposition 2.1. Let Xiea random vector with MEV distribution G, unit Frechet F margins and 
tail dependence function I. If Y has MEV distribution with marginal distributions Fj, j = 1, d, 
and the same copula as X ; then for each A € (0, oo) d , it holds that 

d 

Proof. We first deduce the distribution of \J Fj 1 (Yj). Denoting the copula of X by Cx, we 

i=i 

have, for each u £ [0, 1], 

p ( V m> * «) - & (.*'......*') - o (-^, ^) - 

exp (-/((- log^Ar 1 ,..., (- log^A- 1 )) =u l ( X "'-^ 1) 
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Then 

E ' 



□ 



The next result shows that the max-min coefficient takes into account the tail dependence function 
of all subvectors Xy.^j- of X with indices in Uj^rlj, 7^ T C {1, ...,p}. 

From the Proposition ^. H it holds that, for each 8 / TC {!,■■•, p}, 



' Ev 1 ^) EWii 

e(7, j6 r)^(VV^)|=-^ — V ( 5 ) 



l+ll^AT'^/,),.,^^) 



leading to the following relations of the max-min coefficients with the tail dependence and the 
extremal coefficients. 

Proposition 2.2. //X has MEV distribution then, for each partition T = {7i,...,7 p } of D and 

A G (0, oo) p , it holds that 

iJ(X,A ! l)=e(7 j ,ie{l ) ..,p})- Yl (-l) |T|+1 e(/,-,jeT) (6) 

0#TC{l,...,p} 

and 

fl < x ' 1 - I ) = TT^- E (-i)^' ^^ ■ (7) 

Proof. To obtain the first equality we first apply in the relation 

A V **w = E (-!) |T|+1 V V **w 

j=neij 0#tc{i,..., p } jeneij 

and then the Proposition 2.1 with (|5|). The statement in ([7]) is a consequence of ^ and ([T]). □ 

For the particular case of Ij = {j}, j = 1, ...,p = d, we denote i?(X, A, I) simply by i?(X, A) and 
we have, as a consequence of the above result, that 

l + ^ca..., d} l + i(Ar%(5),... ) AJ 1 <y„(5)) 

and 

e X \ ' , ,s\S\+l £ Xs 



fl (x,i) = -^-- (-i) |s|+ V + 

1 + £X 0^C{T,...,d} 1 + £X ^ 
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These two relations extend the result of the Proposition 1 in Naveau et al. (2009) and equation 
(14) in Cooley et al. (2006), where, for d = 2, 1 - A 2 = Ai = A e (0, 1), we have 



2 ^ ' ' 1 + 1(1^) 2(1 + A)(2-A) 



and R(X,1) - ^=1 



ex + l' 



Our next step is to compare the value of the max-min coefficient i?(X, A, I) with the corresponding 
coefficient in the two boundary cases of independent or totally dependent Xj j: j — l,...,p. 

Proposition 2.3. LetX, X = (Xi, Xj) and X = (X\, Xd) be vectors of unit Frechet random 
variables with MEV distributions such that Xi j , j = 1, ...,p, are independent, X/. , j = 1, ...,p, are 
totally dependent and. for each j = l,...,p, X/., X/. and X/. are identically distributed. Then, 
for each A G (0,oo) d and partition I of D, it holds that 

(a) R(X,X,1) < R(X,X,1) < R(X,X,I), 

(b) R(X,1,1) = 0, 

p 

(c) R(X,l,l) = ^- p £ (_i)|T| + i^_ . 

i+Vex 0#TC{i,..., P } 1 + zL £x ' 3 

3=1 

Proof. If X has MEV distribution then it is associated (Marshall and Olkin (1983)) and then 
the variables M(Ij), j = l,...,p, are also associated (Esary et al. (1967)). For these associated 
variables it holds that 

p j p \ p 

f\ P(M(Ij) > Xj ) > P f]{M(I,) > Xj } > P(M(Ij) > Xj ) 

3 = 1 = 1 / 3=1 

and 

/\ P (M(I 3 ) < Xj ) > P I f] {M{Ij) < Xj} ] > f[ P (Af (/,) < Xj) . 

By taking M^-) = \J X t and M(/ 3 ) = \f Xi, j = 1, ...,p, we can rewrite the above inequalities 

as 

^ fpWi) > > P (f]{M(Ij) > a^j > P ffllM^) > ^} 

and 

P f H Wi) < > P lf){M(Ij) < Xj}) > P I f){M(Ij) < Xj} 
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This concordance order implies (Shaked and Shanthikumar (2007)) that 
and 



\3 = 1 / / \ \j =1 / / \ \J = 1 



for all non-decreasing functions /. The result in (a) then follows by taking / = F and replacing 
X; by p, i £ Ij, for each j = 1, ...,p. 

The equalities in (b) and (c) follows from ((7|) and, in particular for (b) , we recall that if Xj, , 
j = 1, ...,p, are totally dependent vectors then the copula of X is also the copula of the minimum 
(Nelsen, 2006). □ 

For the particular case of Ij — {j}, j = 1, d, the equality in (c) leads to 

fl(xi) = — - v (-i)i s i+ i Iff! = -J— - V(-D fc+i a) — — = — - 

[ ' d+l ^ { ' 1 + 151 d+l J + l d+l 

0^SC{1,. ..,<£} 1 1 fe=l 

which extends the result for the case of d — 2, where i?(X, 1) = i. 



3 Examples and an application 

In order to illustrate the previous results, we consider two families of Multivariate Extreme Value 
distributions and we compute the expressions for i?(X, A, I), which can be easily implemented. 

Example 3.1. //X is the ME V marginal distribution of the Multivariate Maxima of Moving Max- 

oc oo d 

ima processes considered in Smith and Wissman (1996) then l(x\, Xd) = /J /J y x j a l,k,j- 

l—l k— — OO j — 1 

For this tail dependence function we obtain 

«<X,A.X>= Z ( - 1,m+ ' 1 + < - 1)P 



oo oo oo oo p 

MrgiW} 1 + £ £ V V V 1 «W 1 + E E V V K 1 ^ 

1 = 1 k=-oo t£Tjel t J=l fe=-oo t=i jei t 

and, in particular, 

/?rx ai V (-i) |s|+1 i + □ 

«(X,AJ- 2^ 55 55 5 - U 

^ £{ x,..., d}1+ ^ £ i+E E Vv 1 ^ 

d 

Example 3.2. For i/ie Symmetric Logistic model l(x±, ...,Xd) = I I we /iai>e 



\t£T 
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i?(X,A) = £ 



and 



l + (-l) c 



vie 5 / \i=i 



i?(x, i) = J2(-D k+1 (2) t^f - + (- 1 ) d )lT^' 
fc=i 



Several parametric and non-parametric estimators for the tail dependence function are available in 
the literature (Beirlant et al. (2004), Schmidt and Stadtmuller (2006), Krajina (2010)) which can 
be applied to the terms in ©. The comparison of estimation procedures is out of our purposes in 
this paper and we simply remark that the definition of the max-min dependence coefficient suggests 
a non-parametric estimator based on sample means. 

Let X( fe ) = (X[ k , -Xj ), k = 1, ...,n, be a sequence of independent copies of X and Fj the 

(k) 

empirical distribution provided by Xj , k = 1, ...,n, j = 1, a. 
A natural estimator for i?(X, A, I) is 

-, n I p p \ 

£(x,a,z) = 1 53 \/ V $V< W ) - A V > 

fe=i \j : - / j=nei 3 J 

and, in particular for p = d, 

-. n I d d \ 

fe=i y=i i=i / 

If we denote MkiljY 3 = \J F* 3 (A 4 (fc) ) then we can write 

A(X,A,i) = if; \/ M fc (/^ - £ H) |T|+1 ^V 

n k=lj = l 0#TC{l,...,p} fc=lj£T 

The strong consistency of the terms of this sum is stated in the proof of the Proposition 3.8 of 
Ferreira and Ferreira (2012b) and the asymptotic normality can be deduced from the Theorem 6 
in Fermanian et al. (2004). 

As an application of this estimation procedure we consider for X some financial stock markets 
grouped in the tree big world markets I\ = Europe, I2 =USA and ^3 =Far East, as considered in 
Ferreira and Ferreira (2012b). The data are monthly maximums of the negative log-returns of the 
closing values of the stock market indexes CAC 40 (France), FTSE100 (UK), SMI (Swiss), XDAX 
(German), Dow Jones (USA), Nasdaq (USA), SP500 (USA), HSI (China) and Nikkei (Japan), 
from January 1993 to March 2004. 

n 

In the table 1 we present the estimates corresponding to M(A) = i ^3 \/ Fi(Xj) which we need 

k=l i£A 

to compute i?(X, 1,1) for I = {Ii, I2, I3}. 



We obtain i?(X,l,X) = 0,321, R((K h , X/ 2 ), 1, {h, I 2 }) = 0,172, i?((X /x , X/ 3 ), 1, {I u / 3 }) = 
0, 222 and A((X /a , X/ 8 ), 1, {h, I 3 }) = 0, 247. 
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A 


h 


h 


h 




huh 


huh 


huhu h 


M(A) 


0,691736695 


0,614005602 


0,625910364 


0,738655462 


0,770028011 


0,743557423 


0,80070028 



Table 1: ii = Europe, 7 2 =USA and I 3 =Far East 
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